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1. Introduction 

It has been a long held hope [TJ [2] that the large iV limit of QCD is 
simpler than QCD with N = 3 but an analytical solution in d = 4 is still to 
come. A simplification specific to d = 2 resulted in an analytical solution 
to the large N limit of QCD in d = 2 [3]. 

Relatively recently, one has been able to show that continuum reduc- 
tion [1] holds in d = 4 and this has been used to make progress in a numer- 
ical solution of large N QCD. In order to understand continuum reduction, 
it will be useful to understand generalized Eguchi-Kawai reduction [5j . This 
is best motivated by studying the large N limit of two-dimensional U(N) 
lattice gauge theory following Gross and Witten [6]. 

The lectures start off with a pedagogical explanation of the large N 
gauge theories on a two dimensional lattice. We will not include fermions 
anticipating that they do not play a dynamical role in N — > oo limit. After 
this, the generalized Eguchi-Kawai reduction will be explained and it is a 
simple extension of the original argument by Eguchi and Kawai [5] . We will 
then show that reduction to a single site works only in d = 2. Although 
reduction to a single site does not work in d > 2, reduction to a finite 
physical volume will work and this is explained in section [5j 
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The infinite volume limit at finite N is replaced by a infinite N limit 
at finite volume. Spontaneous chiral symmetry breaking can be realized at 
finite volume. Details pertaining to fermions in the large N limit form the 
last part of the lectures. 

2. Large N gauge theories in two dimensions 

Consider U(N) gauge theory on an infinite d dimensional lattice defined 
by the partition function 

Z = I II II dU ^ eS -> S = bNY J Tr (U p + E/J) , (1) 

X fj, p 

where S is the Wilson action, x labels a site, \x labels a direction and p 
denotes a plaquette. U p is the parallel transporter around the plaquette. 
The inverse 't Hooft coupling is denoted by b = p7\f- Under a local gauge 
transformation, 

U X)li -» g x U Xtfl gl + ~ (2) 

and the action is invariant. 

We follow Gross and Witten |6j and gauge fix the two dimensional theory 
by going to the At = gaug43. This corresponds to setting U X) i = 1 for all 
x. In this gauge, 

E TrC/ P = E Tr (^ 2C/ i + i, 2 )- ( 3 ) 

p X 

We still have a remnant gauge symmetry which corresponds to gauge trans- 
formations that are independent of x%, namely g X2 , which can be fixed by 
setting 

U X:2 \ Xl =-oo = !Vx 2 . (4) 
If we make the change of variables, 

U x+ i, 2 = U X U X , 2 , (5) 
then the partition function becomes 

Z=[ l[duJ N ^H U ^) =U[ duj™^), (6) 



This can be done on a infinite lattice but not on a finite lattice. We will address this 
point in section [4] 
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and it factorizes with plaquettes being independently and identically dis- 
tributed. The only remaining symmetry is a global U (N) symmetry, yet © 
is invariant under 

U x - V X U X V% (7) 

for any V x . This is special to Yang-Mills action in two dimensions and is 
not the original gauge symmetry defined in ([2]). Note that ([7]) along with 
([5]) implies ([2]) only if V x = g for all x and this is the remaining global 
symmetry. If we had included fermions, they would have coupled to U x 2 in 
our gauge and we will not have ([7j) as a symmetry. 

The factorization of the partition function enables us to compute expec- 
tation values over individual factors and use that result to get any general 
expectation value. In particular, 

fdUU ije bNT < u+ut ) _ JdUdV(VUV%e bNT < u+ut ) 



JdUe bNTl ( v+u ') JdUe bNTr ( u+ui ) 
jj(TiU)5 i:i = w(b,N)8 tr (8) 



We have used (|7|) in the second equality; 

WVijVtt = ±8 u 6 jk (9) 

in the third equality and w(b, N) is the expectation value of a single pla- 
quette. 

If we define 

z(b,N) = J dUe bNT < u+ut ) (10) 
as the single plaquette partition function, then 

w(b,N) = ^±]nz(b,N). (11) 

The integral in (|10p can be performed [7] and the result is 

z{b,N) = detM; M itj = h-j{2Nb)\ i,j = l,---,N. (12) 

Consider a rectangular L x T Wilson loop with corners at x, x + Tl, 
x + L2 and x + Tl + L2. The parallel transporter around this loop in our 
gauge is 

W X (L, T) = U X+T ^ 2 U x+T i +l ^ 2 ■ ■ ■ ^ r ;r+ T'i+(L-l)2,2 

U x+(L-l)2,2 U x+(L-2)2,2 ' ' ' U xf ^ 
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It follows from ([5]) that 

^ + Ti+(L-l)2,2^J + ( L _i)2, 2 = U x+(T-l)l+(L-l)2 U x+(T-2)l+(L-l)2 

U x+l+(L~l)2 U x+(L-l)2- 

Using dH) and averaging over all the U x variables appearing in the above 
equation, we find that 

l(Tr W X (L, T)> = [e(b, N} T I (Tr W X (L, T - 1)) (15) 

Repeating the above steps L times we arrive at 

±TrW x (L,T) = [e(b,N)] LT , (16) 

since W x (L,0) = 1. The above equation says that the area law is exact in 
two dimensional Yang-Mills theory for all values of N. 

The continuum limit at a fixed N is obtained by taking b —* oo. If 
we first take iV — > oo at a fixed 6 and then take 6 — > oo, we obtain the 
continuum limit of the large N gauge theory a la 't Hooft. In the large 
iV limit, it is instructive to solve for z(b,N) using the method of steepest 
descent resulting in the stationary condition 



2b sin «j 

N 



1 \ - 

V ^— < 



OLi 



(17) 



for the eigenvalues on of the U(N) matrix, U, appearing in (jlOp . Due to 
the symmetry ([7]) of the single site partition function, all expectation values 
will be only functions of the eigenvalues of U and can be evaluated by 
substituting for on the values that solve (|17p. Since the single site partition 
function is dominated by the stationary point, it follows that expectation 
values factorize in the large ./V limit: 

(F(U)) = F((U)). (18) 

The partition function on the infinite lattice can be further reduced from 
a product of independently and identically distributed plaquettes in (|SJ) to 
a single site partition function as in (|1U|) . Expectation of the Wilson loop 
operator defined in (fT3|) reduces to the folded operator 

W x (L,T) = (TrU LT ), (19) 



and it is clear from (|18|) that we will obtain (|16p . This is the motivation 
behind the idea of Eguchi-Kawai reduction [5] . 
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Before we proceed with a discussion of Eguchi-Kawai reduction, it will 
be useful to finish this section with a property of (fTTJ) . We can replace, on, 
by a continuum function a(x), x £ [0, 1] in the large N limit. Furthermore, 
we can define the density of eigenvalues, 

„ W = % (20) 

Then <JT7D reduces to 

2bsma = p[ dpp(fl) cot (21) 

an equation for p where P refers to the principal part of the integral. This 
equation is solved in [6J and the result is 



p ( a ) = I v cos f\/^- sin2 f iib > \ and H < 2sin x yjk ( 22 ) 

[ 27 (1 + 26 cos a) if b < | and |a| < 7r 

The lattice theory undergoes a phase transition at b = \. The continuum 
theory does not exhibit this phase transition. But the lattice strong coupling 
limit and the weak coupling limit are separated by this phase transition. In 
order to obtain the correct continuum limit of the large N theory, we need 
to keep b > ^ and p(a) has a finite region of support around a = that 
does not extend up to a = ±7T. 



3. Generalized Eguchi-Kawai reduction 

The discussion of the large N limit of Yang-Mills theories in section [2] 
suggests: 

1. Factorization of observables; 

2. Domination of the path integral by a single classical configuration. 

Witten argues for the above two points in [8]. Consider, for example, an 
observable that we encountered in section [21 namely, (Tr U x Tr U y ). This 
quantity, in perturbation theory, has two pieces, connected and discon- 
nected. (TrU x ) is of order iV 2 since there are iV 2 gluon species that can 
run around the loop. Therefore, the disconnected piece is of order iV 4 . The 
connected piece has only one loop with two insertions, one for Tr U x and 
another for TvU y , and therefore it is of order iV 2 . Therefore, 



(TrU x TrU y ) = (TrU x ){TrU y ) 



(23) 



G 



in the large N limit. The same argument would also imply that 

^X>0*- (24) 



< e 2 > - < e > 2 = 0; e 



If fluctuations go to zero in the large N limit, it is not necessary to do a 
path integral since one classical field configuration must dominate as was 
seen using steepest descent in section [2l The above argument of Witten was 
made rigorous in [5] and the following statement is a generalization of the 
Eguchi-Kawai reduction: 

Consider U(N) Yang-Mills gauge theory with Wilson action given by 
(TTJ) on a finite lattice of size L\ x L2 ■ ■ ■ with a fixed lattice coupling 
b = p^\f and periodic boundary conditions in all directions. Also consider 
another theory with only one difference from the previous one: = 00 for 
i = 1, ■ ■ ■ , k < d. Now consider an arbitrary closed Wilson loop operator. 
The operators associated with the same Wilson loop on the finite lattice 
and the lattice with k infinite directions could be different due to possible 
folding. The folding comes from the use of periodic boundary conditions 
on the finite lattice. The large N limit will be the same in both cases 
provided the Zn symmetries associated with the Polyakov loops in the fxf, 
i = 1, • • ■ , k directions are not broken on the finite lattice. 

We provide relevant steps for a proof of the above statement by following 
the steps in [5]. Consider a closed Wilson loop that contains the link U X)fl 
once and let us write it is as 

W = Tr U x ^C x>fl (25) 

where C x „ is an open path with more than one link that connects x and 
x + ft and does not contain U x ^. The terms in the Wilson action that 

contain U x u can be written as Tr ( U x U S X u + S x U U X u ) where S x u is the 



sum of the parallel transporters over all the three link paths that connect x 
and x + fi. S Xt ^ does not contain U XjfJi if none of the finite directions are of 
unit lengttH. The group measure is invariant under a small change of the 
form 

U x ,^e leT 'U x ^ (26) 
where T 3 is a group generator and e is a small parameter. Therefore, 



(Tr (T1U X ,,CQ) = ( Tr [T^U^Cl^ 

W Tr( [e" Ti -l] U^St^+S^Vlj [e~^T 3 _ x j ) 



Equation (|30[) remains unaltered if 5^ lM contains U x 
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and to the lowest order in e, 

+bN(Tr (r^U^cQ Tr {tW x ^ - S^E^T'') ) = (28) 
Summing the above equation over all values of j and using the identity, 

TV 2 

= Wfa, (29) 

we get 

(Tr (O*,^) ) + 6(Tr (u^C^U^sQ ) - 6(Tr (c^S^) ) = 0. (30) 

The above equation, referred to as the Schwinger-Dyson equation for Wilson 
loops, relates the expectation value of the original Wilson loop to other 
Wilson loops that correspond to the modification of the original Wilson 
loop by attaching the various plaquette parallel transporters that contain 
U x Even if we start with a Wilson loop where each link occurs only once, 
the above equations will generate Wilson loops where certain links appear 
more that once. In fact, this can already been seen in the second term of 
(|30p . If we repeat the above procedure starting with a Wilson loop where 
U x ,n appears twice, namely, 

W = TrU^C^U^D^ (31) 

where C x ^ and D x „ are two open paths that connect x to x + (1. In this 
case, the Schwinger-Dyson equation will have an additional term of the form 

^(Tr (U X ^CQ Tr (u^Q). (32) 

This is an expectation value of products of Wilson loops. Noting that expec- 
tation value of an closed loop is of order N, we see that the new term is the 
same order in iV as the ones in (|30p and we also note that the new term fac- 
torizes in the large N limit. The coupled set of infinite number of Schwinger- 
Dyson equations obtained in this process will also involve Polyakov loops 
in the finite directions for the following reason: Since it will involve Wilson 
loops of arbitrary size, it will contain loops of the form 



W = TrU x ^U x+L ^Dl^ 



(33) 
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where C x ^ is an open path that connects (x + L^jj) to x + /t and Dj;^ is 
an open path the connects x to (x + (L M + Since Ux+l^^^ = [/^ by 

periodic boundary conditions, this will result in a term of the form 

l(Tr (t^CtJ Tr (u x+l ^dQ ) . (34) 

The difference between the infinite set of coupled equations in the two cases, 
one with k infinite directions and the other being finite in all directions, is 
the presence of additional Polyakov loops in the k finite directions. Polyakov 
loops in the \i direction are not invariant under a global Zjv symmetry in 
that direction where we replace all U x u in a fixed hyperplane perpendicular 

to [i by e l ~U x ^ with < k < N. Since this is a symmetry of the gauge 
action, the Polyakov loops appearing in (|34p will have zero expectation value 
if the Zn symmetry in the \x direction is not spontaneously broken. This 
completes our discussion of the statement concerning generalized Eguchi- 
Kawai reduction. 

4. Reduction to a single site 

The arguments presented in section show that one can reduce the 
large N theory from an infinite lattice down to a single site lattice if the Z/y 
symmetries on a single site lattice are not broken. Since we independently 
showed that this reduction was possible in two dimensions in section [2j it 
follows that the two Zjv symmetries are not broken on a single site lattice in 
two dimensions for all values of b and therefore also in the continuum limit. 
On the other hand, the Z/y symmetries are broken in the weak coupling limit 
in three or more dimensions and we will present the argument following |9j . 

Consider the U(N) Wilson gauge action on a single site d dimensional 
lattice, namely, 

d 

S E K = bN Tr 

The action depends on d U (N) matrices and the gauge transformation is 

U^gU^. (36) 

Note that the eigenvalues of are gauge invariant. We cannot fix a gauge 
such that one of the 17^ = 1 since we are on a finite lattice. The action has 
an additional U d (l) symmetry given by 



(35) 



(37) 
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with — 7r < <j)n < it. The four Polyakov loop operators given by 



(38) 



are gauge invariant but not invariant under (|37|) , If the U d (l) symmetry is 
not broken, then the eigenvalues of all are uniformly distributed on the 
unit circle and = 0. In order to see if this symmetry is spontaneously 
broken in the weak coupling limit, we set 



(39) 



and expand to the quadratic term in the hermitian matrix a^. We fix the 
gauge by setting a± = 0. The group measure is given by 



\\dU, 



fi=2 i>j 



(40) 



where 



$ = sin 2 - [e; - oQ . 



The quadratic piece of the action is 



(41) 



S = -32WV E E a f [P> ij5 ^ ~ fX) a 

i<j fi,u=2 



V ! 



(42) 



where p 13 = p^f . The result of the integration over , ignoring normal- 
ization factors, is 

2-d 



H i>j 



1 1 



W 3 

i>j 



(43) 



Therefore, up to second order in a M , the partition function is 



II IK 



(44) 



All 6 l = has the maximum probability for d > 2 since p 1 - 7 = for all i and 
j implying that the Zn symmetries are broken in the weak coupling limit 
if d > 2. 



10 



5. Continuum reduction 

Consider the continuum theory in a finite torus of size l\ x I2 ■ ■ ■ Id ob- 
tained from the theory on a L\ x L2 • • • periodic lattice at a fixed coupling 
b and taking the limit L\, L2, ■ ■ ■ , Ld, b — > 00 such that 

Li 

li = —=; i = l,2, in a = 2 
\b 
Li 

k = — ; z = 1,2,3 in a = 3 

Zi = 1*0(6); 1 = 1,2,3,4 in d = 4 (45) 
are kept fixed. The lattice spacing, a(b), in d = 4 is given by 

1 /48tt 2 6\i^ 2a^H , 
a(6) = A(~j e_ ~' (46) 



in weak coupling perturbation theory |10| . 

In two dimensions, the variables, 9^ are uniformly distributed in the 
weak coupling limit as can be seen by setting d = 2 in (|44|) . Therefore the 
single site lattice L\ = L2 = 1 will give the same results in the large iV 
limit as any L\ x L2 lattice for all values of the coupling b implying that the 
continuum theory will be independent of l\ and I2. Since the (U(l) in 
the limit N — * 00) symmetries are unbroken for all l\ and I2, it follows that 
two dimensional large N QCD is in the confined phase for all temperatures. 

The lack of a finite temperature phase transition separating the confined 
phase from the deconfined phase is special to d = 2. Numerical analysis has 
strongly established the existence of the deconfining phase transition for 
finite iVmd = 3[m[I2]andd = 4[I3m[IS]. This transition is expected 
to have a large iV limit both in d = 3 [16] and d = 4 [T7j. Therefore, we 
could not have expected single site reduction to work in d > 2. 

Consider a symmetric lattice L\ = L2 = ■ ■ ■ = L (d > 2) and a 
lattice coupling b\ (L) such that no Zn symmetry is broken for b < b\(L) 
but not all Zn symmetries are unbroken for b > b\{L). We know that 
such a coupling exists since all Z^ symmetries will be unbroken for 6 = 
and all will be broken for b = 00. The theory is in the confined phase for 
b < b\{L) and since we expect the continuum theory to have a confined 
phase b\{L) should approach 00 as L — > 00. Fixing, b < b\(L) we can 
consider a L\ x L2 • • • lattice such that Li > L for all i and the arguments 
provided in section [3] shows that there will be no dependence on Li and the 
theory is in the confined phase. Now consider a L x oo d_1 lattice with a 
coupling of b = b\(L). As L is varied we will remain at the phase transition 
point separating the confined phase from the deconfined phase where the 
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Z]\r symmetry is broken in the direction with finite extent L. Therefore, 
l\ = L/b\(L) in d = 3 and l\ = La(b\(L)) in d = 4 should have a finite limit 
as L — > oo and it should be the inverse of the deconfining temperature. This 
has been numerically verified in d = 3 [I] and d = 4 [18] . The order of the 
deconfining phase transition at infinite N can be obtained by a numerical 
computation of the latent heat associated with the transition. The presence 
of a latent heat will result in a jump in the average value of the action 
density at the transition. A non-zero latent heat has been computed in 
d = 4 [19] by numerical methods and numerical studies are currently under 
way in d = 3 [20J . 

We will refer to the confined phase of the continuum theory as the Oc 
phase. There is no dependence on the physical size of the box in this phase. 
There is a transition from the Oc to the lc phase when one of the directions 
has a length less than l\. This is the conventional deconfined phase and there 
is no dependence on the physical size of the box in the (d — 1) directions. 
Let I < Ii be the length of the direction along which the U(l) symmetry is 
broken. Let {d — 2) of the other directions be finite and let one direction 
be finite. Fixing /, we can vary the length of the second finite direction 
and go from a U{1) symmetric phase in that direction to a U(l) broken 
phase. Let > I be the length of the second finite direction such the 
the U(l) symmetry in that direction is unbroken for lengths larger than 
£2(0 and broken for lengths smaller than IzQ,). The existence of ^(0 in the 
continuum theory has been verified by numerical means using the lattice 
theory for d = 3 [21]. The temperature in the deconfined phase is 1 /Z and 
one can remain in the deconfined phase for all temperatures above as 
long as one keeps the extent of the other (d— 1) directions larger than ^(O- 
If we pick one of the other {d — 1) directions to be less than h(l) two of the 
U(\) symmetries are broken and we refer to this as the 2c phase. In this 
manner we can have the continuum theory in a kc phase with < k < d 
where k of the d U(l) symmetries are broken. One can view the (d — l)c 
and the dc phase as the low temperature and the high temperature phase 
of large N QCD in a Bjorken universe [22] . The complete phase diagram 
has not been mapped out in d = 3 or d = 4. 



6. Fermions 

As long as the fermions are in the fundamental representation and we 
only have a finite number of flavors, Nf, fermion loops are suppressed in 
the large N limit compared to gluon loops since NN f « N 2 as N gets 
large |1 , 2j. Physical quantities associated with the fermionic sector of large 
N QCD can be computed using fermionic observables in a gauge background 
generated using the pure gluonic action. Continuum reduction continues to 
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hold while computing physical quantities in the fermionic sector. 

6.1. Chiral condensate 

Chiral symmetry is expected to be broken in d = 2 and d = 4 in the 
confined phase of large N limit of QCD. The theory in d = 2 is in the 
confined phase for any finite torus. The large N degrees of freedom must 
therefore be responsible for spontaneous chiral symmetry breaking in fi- 
nite volume. Consider the lattice model on a single site. We have two 
U(N) matrices, namely, U\ and XJ%. The Wilson action in (pQ) reduces 



to S = 2bNRe Tr 
lattice splits into momentum 



The fermionic operator on the infinite 

slocks with each block being of the form 
Df (Uie tpi , U2e ip2 ;m q ) where — tt < p\,P2 < vr is the momentum of the 
block and m q is the quark mass. The chiral condensate is given by 

X (b,N,m q ) = ^ dpi J* d P2 ±(TrD f {Ux^\U^;m q )), (47) 

where the expectation value is obtained using the single site Wilson gauge 
action. In order to obtain the chiral condensate, we will have to take N —* oo 
limit before we take the m q — ► limit. In the limit of N — > oo, the two U{\) 
symmetries are not broken and (TrDf {U\e ipi , U2e tp2 ; rriqj} is independent 
of pi and P 2- Therefore, the chiral condensate in the massless limit is given 
by 

£= lim Jim ^-(TrD f (U 1 ,U 2 ;m q )). (48) 

m q — >0 N— >oo 1\ 

The fermionic operator Df(U\, u 2 ; 0) on the single site lattice will have 
2iV eigenvalues and it will have iV paired eigenvalues, ±Aj; i = 1, • • • , N 
if the global topology is zero and the fermionic operator obeys chiral sym- 
metry on the lattice. Let us assume that 0<Ai<A2<---<Aat and 
let p(X\, A2, • • • , Xn) be the joint probability of the N eigenvalues after av- 
eraging over U\ and U2 using the single site gauge action as the measure. 
Chiral random matrix theory [23] predicts the joint probability distribution, 
PchRMT^i) Z2, ■ ■ ■ , zn), when TV" is large and where 

zt = Z(b)N\i (49) 

with £(&) being the only adjustable parameter and is the chiral condensate 
at the lattice coupling b. This prediction has been numerically verified [23] 
and the resulting chiral condensate agrees with the known analytical re- 
sult [251126]. 

Turning now to d = 4, we first note following the discussion in section [5] 
that we need to consider a theory on a L 4 lattice at a fixed lattice coupling 
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b such that b < b\{L) and one is in the Zjv symmetric phase. The fermionic 
operators on the L 4 lattice will be such that 

tfj(x + Ljl) = e ip ^(x) (50) 

since the operator in the infinite lattice will split into momentum blocks with 
the momentum translating into the above boundary conditions on fermions. 
We can replace the boundary conditions on fermions by periodic boundary 
conditions and replace the gauge fields by 

U^x) ^ U^x)^ , (51) 

where we have used a gauge transformation to uniformly distribute e iPM over 
L links. Since we are in the U(l) symmetric phase of the large N gauge 
theory, the chiral condensate will not depend on p^ and we can set them 
to zero and compute the chiral condensate. One can use matching with 
chiral random matrix theory predictions to numerically estimate the chiral 
condensate in the large N limit of QCD in d = 4. The eigenvalues, ±Aj, 
i = 1, • • • , 4iVL 4 , will match with the variables, Zi, according to 

Zi = Z{b)NL 4 \i. (52) 

The chiral condensate, £(&), has been numerically extracted and its scaling 
behavior has been studied as a function of b [24] . 

The fermion propagator in the deconfined (lc) phase will depend on the 
boundary conditions, namely (|50p . in the direction where the Zn symme- 
try is broken. Fermions, in this sense, do play a dynamical role in the lc 
phase and anti-periodic boundary conditions with respect to the value of the 
Polyakov loop in that direction will be favored. This has been numerically 
verified in [27]. Since the chiral condensate does not depend on the temper- 
ature in the confined phase, the chiral phase transition in going from 0c to 
lc will be first order. Restoration of chiral symmetry in the lc phase has 
been numerically verified by the presence of a non-zero gap in the fermion 
spectrum for all values of temperature in the deconfined phase [27] . 

6. 2. Meson propagator 

Let 

M(x) = u(x)T S xz {U»)d{z) (53) 

z 

denote a meson at x made out of two different flavors in some spin represen- 
tation given by V. The meson is defined using some gauge field dependent 
smearing operator, S xy {Uu), that commutes with V and transforms as 



S xz (Uf l )=g x S xz {UM 



(54) 
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under a gauge transformation, d(z) — > g z d{z) and it (a;) — ► it(x)gj. The 
meson propagator in momentum space is 



The first equality above assumes that we have translational invariance upon 
averaging over the gauge fields. We have introduced an integral over q in 
the second equality and the integrand does not depend upon q. The tr in 
the second equality indicates a sum over spin and color indices only. We 
extend this to Tr in the third equality where the sum is now over space, spin 
and color indices. The exponential factors on either side of the smeared 
d quark propagator and on either side of the u quark propagator in the 
second equality are viewed as a gauge transformations and results in the 
gauge transformed fields in the third equality. These factors are thought 
of as momenta carried by the quarks: d quark has a momentum equal to 
<7 + 1 and the it quark has a momentum equal to q— |. This corresponds to 
a meson momentum equal to p and q is the momentum around the quark 
loop in the meson propagator. Due to the U(l) gauge invariance in the 
confined phase, we can replace U^e^^ by J7„ in the third equality and the 
integrand does not depend on q. This results in the final equality in the 
above equation. 

The meson momentum can take any value in the range [— ir, tt] since the 
above equation was derived on an infinite lattice. Momenta that are integer 
multiples of ^ are the ones allowed by periodic boundary conditions on 
the L 4 lattice. Momenta that fill in the gaps between the integer multiples 
correspond to boundary conditions of the form given by (|50|) . 

Pion mass and the vector meson mass as a function of the quark mass 
has been numerically studied using the above procedure. This has resulted 



G(p) = Y J eip(x ' y) (M(x)MHy)} 





e 2 m, 




(55) 
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in an estimate of the pion decay constant [28] and the vector meson mass 
in the chiral limit |29j . 

6.3. Numerical details 

We will focus on d = 4 and present the necessary numerical details to 
perform the computations described in the previous subsections. Gauge 
fields are generated using a combination of Cabibbo-Marinari SU{2) heat- 
bath [30j and SU(N) over-relaxation |18j . A description of overlap fermion- 
scan be found in |33 } 134 ^ 135]. The e function appearing in the overlap-Dirac 
operator is best approximated using the 20 th order Zolotarev approxima- 
tion [37j . The action of e on a vector can be performed using the multiple 
mass conjugate gradient algorithm J3SJ[3U]. The low lying eigenvalues of the 
massless overlap-Dirac can be computed using the Ritz algorithm [40J . The 
trace involved in the computation of the meson propagator can be stochas- 
tically estimated using a single random vector. 

7. Other topics and future directions 

Some progress in the understanding of the physical transition from 
strong coupling (hadron resonances) to weak coupling (perturbative QCD) 
has been achieved by the study of Wilson loops |41j . Much more work 
needs to be done in this topic and there is recent progress in connecting the 
transition in Wilson loops to a chiral transition in two-dimensional fermions 
coupled to the four dimensional gauge field [32]. 

A computation of the string tension using folded Wilson loops has been 
performed in d = 3 [33] and d = 4 [33] . Typically, one obtains the string 
tension from Polyakov loop correlators [36] since large Wilson loops have 
a small expectation value due to two reasons: large area and perimeter 
divergence. The results at large N obtained using folded Wilson loops are 
in good agreement with the results obtained at smaller N using Polyakov 
loop correlations and then performing a large ./V extrapolation [46\ I47j . 
Based on our discussion in section [3l we note that folded Wilson loops can 
also be computed in the lc phase as long as the Wilson loop is in the plane 
where the Zn symmetries are not broken and such loops are referred to 
as spatial loops. Numerical computation of the spatial string tension in 
the lc phase of d = 3 show that the string tension grows linearly with the 
temperature [38] . 

Fermions play a dynamical role in the large N limit when they are in 
the adjoint representation. Since the gauge action induced by the fermion 
determinant will tend to cancel the original gauge action, one expects single 
site reduction to hold in d = 4 if the fermions are in the adjoint representa- 
tion. Strong arguments in this direction have been put forward recently in 
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the continuum |49} 150] . Taking the zero volume limit in the continuum is 
not necessarily the same as working on a single site lattice. Therefore, it is 
interesting to numerically study large N gauge theory with fermions in the 
adjoint representation on a single site lattice. Recent progress on this topic 
can be found in |5 H I52j. 
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